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Abstract: We consider a general framework to study holographically the dynamics of 
fundamental quarks in a confining gauge theory. Flavors are introduced by placing a set of 
(coincident) branes and antibranes on a background dual to a confining color theory. The 
spectrum contains an open string tachyon and its condensation describes the U{Nf)i x 
U(Nf)ji — > U(Nf)v symmetry breaking. By studying worldvolume gauge transformations 
of the flavor brane action, we obtain the QCD global anomalies and an IR condition that 
allows to fix the quark condensate in terms of the quark mass. We find the expected Nj 
Goldstone bosons (for m„ = 0), the Gell-Mann-Oakes-Renner relation (for m q small) and 
the r/ mass. Remarkably, the linear confinement behavior for the masses of highly excited 
spin-1 mesons, ~ n is naturally reproduced. 
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1. Introduction 

The original AdS/CFT conjecture |2[ states the existence of a stringy description of a 
conformal, highly symmetric, quantum field theory. The low energy limit on the string side 
(supergravity) corresponds to the strong coupling regime of the gauge theory. A natural 
question that has attracted a lot of attention is whether such a duality can be generalized 
in order to describe the strong coupling regime of more realistic gauge theories. In the last 
few years we have learned that the ideas of holography can indeed be applied to a certain 
extent in that direction. In particular, higher-dimensional duals of (non-supersymmetric) 
QCD-like theories have been built and many qualitative and, surprisingly, some quantita- 
tive features of real-world strong interactions can be described in such frameworks with 
reasonable accuracy. 

One of the most important strong-coupling phenomena of QCD is spontaneous chiral 
symmetry breaking U(Nj)l x U(Nf)ji — > U(Nf)y by the formation of a non-perturbative 
quark condensate (qq) 7^ (Nf denotes the number of quark flavors). The goal of this 
paper is the study of this symmetry breaking using holographic techniques. 

We start by briefly reviewing the results obtained on this topic in the existing literature. 

This question was initially addressed in (3|, [|] (see also || ) , by considering stacks of Nf 
flavor branes in non-supersymmetric backgrounds created by N c color branes. The chiral 
symmetry is a U(1)a isometry of the geometry and is spontaneously broken due to the 
embedding of the flavor branes. Typically, these models are asymptotically supersymmetric 
in the UV. This allows to have good control on quantities like the quark condensate and the 
quark mass. Imposing regularity of the brane configuration, one obtains an IR condition 
which fixes the quark condensate in terms of the quark mass, as expected in QCD. Also, 
by introducing a small quark mass, the Gell-Mann, Oakes, Renner (GOR) relation for the 
pion masses 

2 n m q {qq) 

Jit 

can be obtained. The 0(^f) mass of the rf in this kind of setup was studied in Q. 
The limitation of this approach is that, even for massless quarks, only the abelian chiral 
symmetry is present, while the non-abelian chiral subgroup of the flavor symmetry is absent 
from the beginning, and its breaking cannot therefore be seen. From the field theory 
side, this is due to the existence of a q&q Yukawa coupling (color and flavor branes are 
codimension four), where $ is one of the (massive) scalars living on the brane world- 
volume ||. 

A different approach was followed in [|], ||. By considering Nf D8-D8 pairs in a 
non-supersymmetric background corresponding to N c D4 branes, a U(Nf) x U(Nf) global 
symmetry was introduced. Contrary to the case described earlier, the chiral symmetry is 
now realized on the flavor brane world- volume rather than in the geometry. The breaking 
of the full non-abelian chiral symmetry is due to a smooth recombination of branes and 
antibranes in the IR. By computing the meson spectrum, Nj massless Goldstone bosons 
are found. One of them, the rf, is massless only when N c — > 00 and the authors addressed 
the problem of finding its 0(j^-) mass. The chiral anomaly and associated WZW term 
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are nicely reproduced. Remarkably, quantities (some meson masses, decay constants and 
couplings) computed in this holographic setup match reasonably well the experimental 
values. Moreover, the model also incorporates a very natural picture of chiral symmetry 
restoration at high temperature: once a horizon is formed in the geometry, the branes and 
antibranes can fall into it and they do not need to recombine |J. 

This construction however presents some shortcomings too. First, there is no param- 
eter which can be associated to the quark bare masses. By modifying the flavor brane 
embedding, one can vary the masses of the massive mesons, which can be thought of as a 
modification of the constituent quark masses. However, the pions remain massless, imply- 



ing that the bare quark masses are always vanishing (see ||, [LC], |TT], [L2|]). This is clearly 
not a physical feature of QCD. Related to this problem, there is no parameter that can 
be identified with the quark condensate, even though this is a crucial quantity to describe 
spontaneous chiral symmetry breaking in QCD. Another unphysical feature of the model 
is the absence of a tower of massive mesons with the quantum numbers of the pion. 

Finally, there is an interesting alternative approach, that has been named AdS/QCD 



|13| , 14]. It is a bottom-up construction in the sense that it does not emerge directly from 
some known string theory. Rather, it just assumes the existence of a dual description of 
QCD living in a five-dimensional, asymptotically AdS space. The five-dimensional bulk 
fields needed to dualize the quark bilinears are introduced by hand, and spontaneous chiral 
symmetry breaking is modeled by giving a vev to a bulk scalar by hand. The model 
incorporates the existence of Goldstone bosons, vector meson dominance [14] and the GOR 
relation for the pion masses Jl3j . In its simplest form, the five-dimensional space is taken to 
be just AdS§ with a hard IR cutoff. However, modifications of the background have been 
considered in order to incorporate the effect in the geometry of the possible condensates jl5| 
or linear confinement [p^ ]. A drawback of the model is that, since the condensate is not 
determined from some dynamical computation, there is an extra parameter compared to 
QCD. Also, being ad hoc, the models are less satisfying and it would be nice to understand 
how they are related to some (possibly non-critical) string theory which may help fix the 
background geometry or the potential for the bulk scalar. Some progress along this line 
seems imminent, [17]. 

The goal of this paper is to improve the present holographic description of chiral 
symmetry breaking by proposing a general string theory construction which accounts for 
all relevant phenomena in a natural and simple way. Flavors are introduced via a stack of 
Nj overlapping brane-antibrane pairs. At the same time, we allow for a non-trivial profile 
for the open string tachyon field. This is generally the lightest string mode extending 
between branes and antibranes, and it transforms in the bifundamental representation of 
the U (Nf)i, x U (Nf)n flavor symmetry group supported by the brane-antibrane pairs. This 
complex scalar is therefore the natural candidate to describe chiral symmetry breaking, as 
suggested in {?]] (see also [11, 18, [0|]). A review on open string tachyon physics can be 
found in iSI. 



We will study the world-volume effective action for the brane-antibrane system, in- 
cluding the tachyon both in the Dirac-Born-Infeld and Wess-Zumino terms. To keep the 
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focus on the relevant dynamics, we will consider a simple setup in which only the tachyon 
and the world- volume gauge fields of the branes and antibranes are dynamical. This is the 
minimal setup to include the quark bilinears in the dual theory (up to spin one) and it 



matches the bulk field content of the AdS/QCD models 14 1. We will also not worry 
about the closed string physics. We will keep the discussion very general, without applying 
it to any particular model: we will only need to assume that the background metric is such 
that the color theory is confining (using a general prescription in 

In a few places we will also perform some more explicit computations; in these cases 
we will assume that the space is asymptotically AdS in the UV. This will allow us to apply 
the general formalism developed in this paper to a relevant and interesting example, while 
hopefully clarifying our presentation. 

1.1 Summary of results 

The open string tachyon we introduce transforms in the bifundamental representation of 
the flavor symmetry group U{Nf)i, x U(Nf)n, and couples on the boundary to the quark 
scalar and pseudoscalar bilinears. Up to a normalization we will introduce later on, we 
have in fact: 

„ - i + 7 5 -i-V* ,, u 

T^q — - — q T' <-> q — - — q (1.1) 

By considering brane-antibrane configurations, and including the dynamics of the com- 
plex open string tachyon scalar T, we have been able to provide a general framework to 
describe holographically the chiral dynamics of QCD and related theories, and to reproduce 
many features of QCD: 

• The UV non-normalizable component of the tachyon corresponds to the quark's bare 
(complex) mass matrix m q . For simplicity, in this paper we will consider all quarks 
to have the same (possibly vanishing) mass m q . 

• The UV normalizable mode of the tachyon corresponds, in turn, to turning on an 
expectation value for the quark condensate. 

• We show that in a confining background, even for m q = 0, the tachyon profile cannot 
be trivial, which results in the chiral part U(Nj)a of the flavor symmetry being 
always broken, either explicitly or spontaneously. The construction of this paper, 



therefore, provides a holographic version of the Coleman- Witten theorem [22|. 



The fact that the tachyon diverges in the IR (signalling the IR fusion of branes and 
antibranes) constrains the way it vanishes in the UV: this allows to fix the value of 
the quark condensate (qq) in terms of the quark bare mass m q . 

We can derive formulae for the anomalous divergences of flavor currents, when they 
are coupled to an external source. 

The WZ part of the flavor brane action gives the Adler-Bell-Jackiw U(1)a axial 
anomaly p3|] , and an associated Stuckelberg mechanism gives an O ( -jJ- ) mass to 



the would-be Goldstone boson rf , in accordance with the Veneziano- Witten formula. 
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• When m q = (i.e. the tachyon has no UV non-normalizable mode), we find that 
(in the iV c — > oo limit) the meson spectrum contains Nj massless pseudoscalars, the 
U(Nf)A Goldstone bosons. 

• Studying the spectrum of highly excited spin-1 mesons, we find the expected property 
of linear confinement: m\ ~ n. 

• When we consider an asymptotically AdS space in the UV, we can precisely reproduce 
the GOR relation for the mass of the pions in the case of a small but non-zero quark 
bare mass. 

On the other hand, the main limitations our approach still presents are the following: 

• The details of the tachyon potential that one should insert in the DBI action are not 
rigorously known. Despite this, its gross features suffice for extracting qualitative 
information. 

• A background in which the brane-antibranes are on top of each other (and with no 
Wilson line turned on), would probably require large curvature, as we explain. Then, 
a reliable background metric cannot be found by solving Einstein's equations. 

• We find that the slope of the linear confinement relation is not the same for vector 
and axial mesons. This may be due to the fact that it is not clear whether the 
DBI action we use to include the effects of the dynamics of the open string tachyon 
remains valid in the region where the tachyon and its derivative diverge. 

However we found encouraging that a large amount of expected QCD features can be 
extracted from our construction despite the misgivings above. 



2. The general construction 

The general setup we consider consists of a system of Nf overlapping Dp-Dp flavor brane- 
antibrane pairs in the background associated to N c Dq color branes. In critical, ten- 
dimensional string theory, we will need to impose the condition p + q = 12 and require the 
color and flavor branes to be codimension 6, in order to get the correct WZ couplings. 1 For 



instance, one can think of a D3-D9-D9 intersection as in [18] or a D4-D8-D8 intersection 
similar to the Sakai-Sugimoto model 0]. Notice, though, that the general formalism we 
develop can also be applied to lower dimensional setups provided the WZ couplings are 
the obvious adaptation of those of section 0. Interesting examples of this kind are the 



1 The reason why is very simple to explain. In section 3.2 we will be looking for a coupling on the flavor 
Dp — Dp branes which involves 6-dimensional gauge field forms. This can only come from the Wess-Zumino 
coupling of the p-branes to a magnetic fc-form potential, where fc = p + 1 — 6 = p — 5. Now Dg-branes 
source a C q +i RR potential, which is parallel to the source branes and generically depends on the radial 
direction alone. Imposing 4-dimensional Lorentz invariance of the background ensures that the Hodge dual 
of C q +i is a magnetic fe-form, with k = 10 — (q + 2) — 1 = 7 — q. Imposing k = k finally gives the condition 
p + q = 12. 
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AdS/QCD models or non-critical string constructions (jV = 1 theories with space-time 
filling brane-antibrane flavors have been discussed in |24], |19L [25|] and M = in p!o|). 

In the spectrum of stacks of overlapping brane-antibrane pairs there is a tachyonic, 
complex scalar, open string mode, which transforms in the bifundamental representation 
of the U(Nf)z x U(Nf)n flavor symmetry group supported by the brane-antibrane pairs. 
As we argued in the introduction, we can in general associate chiral symmetry breaking 
to a vacuum expectation value for this field. Near the UV region (in our coordinates this 
is around z = 0), the tachyon's vev vanishes and the full U{Nf)i x U(Nf)n is present, as 
in the lagrangian of massless QCD. As we approach the IR, we prove in section [T^ that if 
the theory is confining the tachyon must acquire a non-trivial vev, breaking the symmetry 
to its diagonal U(Nf)y subgroup, (see also section |2.2p . Concretely, we show that the 
vev reaches infinity at a finite distance in the bulk, a point that we denote by zir. In a 
sense, this process can be thought of as brane-antibrane recombination. In the minimal 
construction that we consider in the rest of this paper, zir coincides with the end of space. 
Figure [j] presents a general profile for the tachyon in the setup we just described. 

Z 



UV 




z = o z = z IR 

Figure 1: A scheme of the setup. Generically, the tachyon vanishes in the UV where one has a 
stack of brane-antibrane pairs. In the IR, it diverges and the brane-antibrane pairs end smoothly. 
Here r stands for the modulus of the tachyon field. 



For the present discussion, it will be very important to consider the completion of 
the world-volume action of the brane-antibrane system to include also the couplings to 
the open string tachyon. As usual, this action is composed of a Dirac-Born-Infeld and a 
Wess-Zumino term. In the rest of this section and in section [|, we study the physics of the 
DBI piece, which yields the vacuum of the configuration and the meson mass spectrum. 
The WZ part will be studied in section ||. From it, we derive the holographic description 
of the parity and charge conjugation symmetries, the global anomalies, and a holographic 
Coleman- Witten theorem. 



2.1 The degrees of freedom: a minimal setup 



As proposed in [26, p7| 1, the semiclassical action for a Dp-Dp system which includes the 
physics of the open string tachyon is a generalization of the usual Dirac-Born-Infeld action 
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for D-branes. 

We start from the simpler configuration of a single brane-antibrane pair. In this case 
the proposed DBI action reads: 

S = - J d p+1 x e~* V(t 2 , Y[ - x) (y~ det A L + y 7 - det A fl ) (2.1) 

where (we will set 27ra' = 1 from now on, and only reinsert a' when needed): 

A {i)M N = P[9 + B]mn + F$ N + d M Y^d N Y^ + ±-(D M T)*(D N T) + ^(D N T)*(D M T) 
f mn = 9 mA [ S - d N A® , D M T = (d M + iA L M - iA^)T (2.2) 

and V(t 2 ,Y[—Y^ : x) is the tachyon potential. The complex tachyon is denoted by T = re te , 
the indices i = L,R denote the brane or antibrane, the Yl\ are the transverse scalars, and 
A® is the world-volume gauge field. Here and in the following, capital latin characters, 
M, N, denote world- volume directions, which include in particular the Minkowski direc- 
tions, which we will denote by greek letters fj,, u, and the holographic radius which will be 
referred to as z. 

For simplicity, we will not consider a background i?-field. We will also ignore the 
transverse scalars living on the flavor branes. These modes, which are generally present 
in a critical string setup, do not have any obvious QCD interpretation and would not 
play any relevant role in the present discussion. Thus, we only consider the degrees of 
freedom coming from the open string tachyon and the world-volume gauge fields along the 
Minkowski or the holographic directions Aji , Az ■ This field content is enough to include 
the basic QCD quark bilinears up to spin 1 and coincides with the one of the AdS/QCD 
models. 

Higher spin operators correspond to stringy excitations and cannot be addressed in 
this formalism. They could be introduced by hand as in [pg| . 

Since we are ignoring transverse scalars, for the tachyon potential we take the expres- 
sion for overlapping brane- antibranes 2 derived in boundary string field theory [29], [ 30 1 3 : 

V(T 2 )=T p e- T ' 2 (2.3) 

One should keep in mind that this expression is not rigorously justified since there can be 
non-trivial field redefinitions between the two formalisms. Moreover the boundary string 
field theory computation is performed in flat space and it is not trivial at all that the result 
should apply unchanged to curved backgrounds. Remarkably, in section [D] we will find 
the nice feature of linear confinement using this gaussian behavior for the potential. Most 
of the rest of the discussion is insensitive to the precise expression of the tachyon potential. 



2 If one considered more complicated setups with non-trivial brane profiles, the potential should depend 
also on the brane-antibrane distance. For a discussion in a related setup, see 

^Unfortunately, there are not unified conventions in the literature about the definition of T. The un- 
ambiguous statement is that the mass squared of the open string tachyon in flat space is = — — tj. 



Equation (2.2.) is thus consistent with ( |2.l| ) and (2.2) since, then, for small r the action is proportional to 
S oc J (|((9t) 2 — 5 2a 7 '?" 2 ), wnere we have considered flat space and reinserted a' . 
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After presenting the simpler abelian case, we can now go to considering a stack of 
Nf > 1 brane-antibrane pairs, which is what we are eventually interested in. The natural 
non-abelian generalization of ( |2.lD is: 



S = - J d p+1 x SymTr e" *V (TT f )( V - det A L + V" det A R ) (2.4) 



where SymTr is a symmetric trace as defined in [31]. The analogous proposal for the 



action of a stack of non-BPS D-branes was made in [32] and can be related to systems of 



branes-antibranes along the lines of [27]. Although the symmetric trace prescription for 
the non-abelian action is known to be corrected at higher orders in a' , this will not be 
relevant for our discussion. All that is important is that the potential appears in a single 
trace, and this is guaranteed by the large N c limit, [p^] . 

The gauge group supported by the brane-antibrane system is U(Nf)L x U(Nf)n. The 
gauge fields Al,r transform in the adjoint representation of U{Nf)L and U{Nf)n respec- 
tively, while the tachyon T transforms in the bifundamental of U{Nf)i x U(Nf)n, (that 
is, in the antifundamental of U(Nj)l and in the fundamental of U(Nf)n). Its hermitian 
conjugate transforms in the opposite way. For more details on our choice of conventions, 
we refer the reader to appendix |A[ 

2.2 The vev for the tachyon field 

In this section we want to determine the vacuum configuration for a stack of brane- 
antibranes over a confining vacuum. It is trivial to show that the equations of motion allow 
for the gauge fields to be consistently set to zero, (Al) = (Ar) = 0. As we mentioned in the 
introduction, and we will explicitly show in section |3.2] , anomaly considerations regarding 
the WZ part of the action of a stack of brane-antibranes over a confining vacuum require 
that (t)\ zir — > oo. This necessitates that we allow for a non-trivial profile of the tachyon 
in the vacuum configuration. 

Four-dimensional Poincare invariance determines the general form of the bulk five- 
dimensional metric. In the mostly plus signature this reads: 

dsl = g xx (z)r]^ u dx^dx u + g zz (z)dz 2 (2.5) 



The non-abelian action fl2.4| ) is a very complicated object, where different components 
of the fields mix because of non-trivial commutation relations. We start from the simpler 
Nf = 1 case, where these complications do not arise. From (|2.l| ), it is consistent to set the 
tachyon phase to a constant. 

We are thus left with the world-volume action for the modulus of the tachyon r, which, 
upon reducing to five dimensions reads: 



S = -2 J d 4 xdze~^ff V(T 2 )g 2 xx] Jg zz + \{d z rf . (2.6) 
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The prefactor e ^ e ff(z) is defined as the dilaton exponential times whatever comes from 
integrating the transverse, spectator, p — 4 dimensions. Schematically: 



e-^ff(z) V(r\z)) g 2 xx s/g zz + {d z rf = 



(2.7) 



= [ d p - A y e~+ V{t 2 , Y[ - x) ( y/- det A L + V" det A R ) 

where A/" p_3 is the internal part of the flavor branes world-volume. 4 
The Euler-Lagrange equation for r obtained from ( |2.6| ) reads: 

+ 1^^) + ^<M|L£_^l _ 4T + 2T ( * + (4T)2) _ 

t 9zz9i x e ^tt \ g 2 x e <h/f 2g zz ) V 2 J 

(2.8) 

This is a second order differential equation which therefore has two integration constants. 
Since r is dual to the quark bilinear, we know on general AdS/CFT grounds that, by 
looking at the UV behavior of r, these two constants can be related to the quark bare mass 
and the quark condensate. We now assume that the space is asymptotically AdS, that is 
9xx{ z ) — 9zz( z ) — R-AdSi ^ + • • • an d e~^ e ff ~ const = e~^° for small z. The tachyon field 
is dual to the quark bilinear which has dimension 3. Thus, the usual relation between the 
mass of a five-dimensional scalar field and the conformal dimension of its dual operator, 
A(A — 4) = m 2 !^^' nxes rn T^AdS = ~ ^> wn i cn m turn implies that, if = — ^7 = — 7T, 
the radius of AdS should be small: R\ dS = 6a' = j. This leads to the following UV 
behavior (in order to simplify notation the vacuum expectation value (r) will be denoted 
just as r from now on): 

3 3 / 4J 1 \ 2 

Rl dS T m n = R% dS y-^r^) r = m q (z + ...)+a(z 3 + ...) (small z) . (2.9) 
We have written this expression in terms of the canonically normalized r can , such that the 

3 

boundary coupling is R\ dS J d xr can z~ l qq. The parameter a can be related to the quark 
condensate (qq), see appendix || for details. The UV behavior (|2.9[) is the same as the one 
for the scalar breaking the chiral symmetry in the AdS/QCD models [|^, 14]. 



We now want to prove that the IR condition t\ Zir — * oo gives a condition on t(z), 
fixing a, and therefore the condensate, in terms of the mass m q . 

A general result of the gauge-gravity correspondence states that a sufficient condition 
for confinement is that at some z div , g zz (z div ) -> oo while g xx {z div ) / and d z g xx (z div ) < 



1 21]. We will assume this is the case for our background and in particular, after possibly 



redefining z, we fix the divergence to be a single pole g zz = b(z d i v — z) 1 near z = z d i v . 



4 The reader might be puzzled by the approximation involved in discarding all dependence on the in- 
ternal world-volume scalars. It should be noted, though, that our ansatz for the five-dimensional metric 



(2.5) is completely general, being constrained only by four- dimensional Poincare invariance. The vacuum 
configuration of the internal scalars will only affect the explicit form of the metric coefficients and the 
radial dependence of the tachyon modulus and dilaton. Since our arguments are completely independent of 
the particular form of the confining, reduced, five-dimensional background, all the following results apply 
regardless what the critical string setup is. 
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This is the behavior, for instance, in Witten's dual of Yang-Mills Inspecting the 
leading terms for large r of fl2.8| ), we see that the tachyon can only diverge at zjr = Zd 



IV 



is a positive 

ZIR 



and there it diverges as r oc [ztr — z) a , where a = —irb 

number since the derivative inside the bracket is typically negative 

We warn the reader that this result should be taken with a grain of salt because it relies 
on the behavior of fl2.g| ) when r, f — * oo, which lies outside the regime of validity of the DBI 
action (p.6|). Nonetheless, the general lesson one should learn from this computation is that 
the t\z ir — ► oo consistency condition constrains the IR behavior of the tachyon. Physically, 
this fixes the quark condensate in terms of the quark mass 5 , in the spirit of |3], Q . 

At this point we can consider the more general non-abelian case Nf > 1. Throughout 
the paper, we will consider for simplicity that all Nf flavor fields have the same mass m q . 
The symmetrized trace in the action (|2.4| ) might seem to make it a lot more complicated 
to find a vacuum configuration for the tachyon matrix T. Fortunately we can show that 
there is a perturbatively stable configuration given by Nf identical copies of the solution 
for the abelian case we described above 

(T) = t{z) I (2.10) 

In fact, we can apply here a result that was derived in 0] for a different setup. The 
reasoning in showing that the vacuum configuration where vevs are proportional to 
identity matrices, like fl2.10|) in our case, is a stable (up to quadratic order) solution of the 
non-abelian DBI action is very general. It only relies on having a single trace in front of 
the action and the fact that the difference between abelian and non-abelian DBI actions 
consists entirely of terms involving commutators. The argument does not prove that this 
is the global minimum. There could in principle be a different configuration with lower 
energy, but this is an unlikely possibility and we will assume that ( 2.1C| ) describes the 
vacuum of the theory. 

Once we have established that (T) is proportional to the identity matrix, the equation 
of motion for r(z) is the same as ( [2.8D as well as the UV and IR conditions, since for 
the configuration ( p. 10 ) and up to the order we need here, the non-abelian action (|2.4| ) is 



actually the sum of Nf identical copies of the abelian action (2.1) 



Looking at how the fields transform under flavor symmetry transformations ( |A.6| ), it 



is clear that the vacuum ( 2.10 ) spontaneously breaks U(Nf)i x U(Nf)ji into its diagonal 



subgroup U(Nf)v, where the vectorial gauge transformation is Vl = Vr. 



It is natural that for a fixed m q , the conditions stated above fix uniquely a. However, within the present 



iA 



very general formalism, this uniqueness cannot be proved. We will come back to this in section 

6 From the field theory point of view the nature of the breaking depends on the asymptotic UV behavior 
of the tachyon. When m q ^ 0, the vacuum configuration of the tachyon is not normalizable and we are 
adding a mass term for the quarks to the lagrangian: the breaking is explicit. If instead m q = 0, the 
tachyon is normalizable, and we are turning on an expectation value for the quark bilinear in the original 
theory, which makes the symmetry breaking spontaneous. From the gravity point of view, instead, both for 
vanishing or non-zero m q , the breaking is spontaneous since it is determined by a non-trivial expectation 
value. 
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When niq = 0, this spontaneous breaking implies the existence of iV? Goldstone bosons. 



We will explicitly find them in the spectrum in section pL3| (see also section 3.5 for an indirect 
argument). We conclude this section by noting that the fact that the symmetry is broken 
down to U(Nf)v, and not further, is also related to a theorem by Vafa and Witten [34] 
that states that vectorial symmetries cannot be spontaneously broken. 

3. The D-brane Wess-Zumino sector 

We now consider the WZ coupling of the Nj Dp -Dp branes to the background RR fields, 
including the tachyon dependence. This part of the action is given by (the following 



expression was proposed in |35f| and proved in |3j], 37] using boundary string field theory): 

S wz = T p J C A Str exp [i2vra'^] (3.1) 

where is the world- volume of the Dp -Dp branes, C is a formal sum of the RR potentials 
C = J2 n (—i) E ~^ C n , and T is the curvature of a superconnection A. In terms of the 
tachyon field matrix T and the gauge fields Al and Aft living respectively on the branes 
and antibranes, they are (from now on, we will set 2na' = 1 again and use the notation of 

In appendix |c] we review the relevant definitions and properties of this supermatrix for- 
malism (in particular, notice the definition of the supermatrix product). 
The superconnection is defined as: 

T = dA - iA A A (3.3) 

and satisfies the Bianchi identity: 

dT - iA A T + iT A A = (3.4) 



Using this identity and the cyclic property of the supertrace ( |C.3| ) it is immediate to 
check that Str e %:F is a closed form and therefore that ([O]) is invariant under RR gauge 
transformations. At least locally, there exists then a form £1 such that: 

dn = Stre^ (3.5) 

The case of interest here is the WZ coupling on the Dp -Dp flavor branes in the back- 
ground of N c Dg-color branes. As explained at the beginning of section B we have p+q = 12 
and the electric RR potential sourced by the color branes is then Cq+i = Ci3_ p , while its 
magnetic dual is C p _5. Formula ([O]) reduces to: 

Swz = iT p I C p _ 5 A Str e^| 6 _ form = i{-f ' T p I F p „ 4 A ft| 5 -form (3.6) 
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When the branes worldvolume has boundaries, the two expressions in (|3.q ) differ by 

boundary terms, and are therefore not equivalent. In this case, which is the relevant one 
for this paper, it is argued in |38| (for intersecting branes) that the correct form for the 

action is the last one in fl3.6| ). Since 1 5 form is defined only up to an exact form, there is 

then a possible ambiguity in the definition of the WZ action Q3,6| ), We will assume that 
the correct form for f2|5_f orm respects the discrete symmetries of Str e ljr |6-form> which we 
will discuss in the following subsection |3.1| . 

Notice that the potential sourced by the Dq-branes has indices parallel to the branes 
and generically depends on the radial direction. Moreover, 4-dimensional Lorentz invari- 
ance of the background (plus possibly a radial redefinition) prevents the metric from having 
non-zero off-diagonal terms in the directions (0, 1, 2, 3, z). These two facts ensure that F p _4 
does not have legs along the (0, 1, 2, 3, z) directions. Thus, the only components of the five- 
form f2|5_form = ^5 which we are interested in are those along X .. X . X , X and z (we 
will call M§ this five-dimensional space). Similarly to section ||, we assume at this point 
that does not depend on the directions along which F p _4 lies, as expected on physical 
grounds. Then we can integrate out those directions in ( |3.6| ) to obtain the 5-dimensional 
effective Chern-Simons action: 

iN c 



Scs 



p+i 



nMi 



•p-4 



Ms 



(2vr) 5 a 



(3.7) 



M 5 



where we have used the quantization condition f = N C T\2- P and the explicit 



value of the D-brane tensions. 



3.1 Discrete symmetries (P,C) 

The six-form Str e^|6„f orm appearing in ( |3,€| ) is invariant under two discrete symmetries, 
P and C, that we will present in this subsection. Thus, an invariant under P, C exists. 



It follows then from (3/7) that the CS action is invariant under the same symmetries, which 
acquire a fundamental physical interpretation in the dual field theory. 

Parity 

The six-form Str e l ^\s_{ ovm is left invariant by the change P = P\ ■ P2, where: 

P 1 :A L ^A R , T^T ] and P 2 : (a?i, x 2 , x 3 ) (-xi, -x 2 , -x 3 ) (3.8) 
Notice that Pi also transforms DT <-> DT\ Fl <-> Fr. Indeed we can easily check that: 

Pi 



Str e 



Str 



cxp 



iF L - T^T DT^ 
DT iF R - TT^ 



Str 



cxp 



iF R - TT^ DT 
DT* iF L - T^T 



-Str 



exp 



iF L - T^T DT^ 
DT iF R - TT^ 



. (3.9) 



For the last equality we have just swapped the ordering of the blocks inside the matrix, 
picking up a minus sign due to the definition of the supertrace (p.2[). The overall minus 
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sign in (3J)) is compensated by the minus sign acquired by the form dx\ A dx2 A dx% under 



the action of P2, and the CS action (3.7) is invariant under P. 

We conclude that four-dimensional parity is equivalent to charge conjugation for the 
D-branes (interchanging D <-> D). 

Charge conjugation 

The other discrete symmetry of Str e 4jP |6_f orm , and consequently of the CS action 

t 



is: 



C : 



A, 



-a* l , t — > r* 



(3.10) 



where the superscript t denotes matrix transposition. Notice that under ( |3.10| ), the field 
strengths and tachyon covariant derivatives transform as Fl — ► — F R , Fr — > —F L , DT — > 
DT 4 , DTt — > (DT^Y . The supertrace transforms as: 



Str e 



Str 



exp 



Str 



exp 



-Str 



-iF R -(TT^y (DT*) 

DT t -if£ - (T^Tf 

iF L - TtT -iDTt 



-Str 



/ - (T T T)* 

6XP 1 (DTt)' -if* -(TTt)* 



exp 



-iDT 



iF R - TTt 



pt" 



-Str 



exp 



-if L - T^T -iDT^ 



-iDT 



-iF R - TTt 



(3.11) 



For the identity in the second line we have again swapped the ordering of the blocks 
picking up a minus sign. In the equality of the following line we used the definition of 

and the fact that the 



pseudotranspose (C.4). Finally in the last one, we used property flCT 
supertrace does not change under pseudotransposition. The matrix of the last expression 
is not the same as the untransformed one of the first line of (3.11). Therefore, fl3.10| ) does 
not leave the supertrace invariant. However, we are only interested in the 6-form within 
the full expression of the supertrace. Noticing that all two- forms (Fl, Fr, DT ADT') inside 
the matrix of the last line have picked up a minus sign compared to the initial expression, 
we find that the 6-form (in fact, all the 4k + 2-forms) is left invariant. 

If one wanted to maintain invariant the 4fc-forms of the supertrace (which can be 
related to gauge theories in 4A; — 2 dimensions in the same way as we relate the 6-form to 
a four-dimensional field theory), then the transformation would be Al — > —At , Ar — > 



in Ak 



(T^y , Tt — r T t . This is in agreement with the fact that charge conjugation 
2 dimensions does not change chirality. 



Table 1 summarizes the transformation properties of the different degrees of freedom. 
In particular, notice that the vacuum condensate ( |2.10 ) transforms as ++ as expected for 
the (qq) bilinear. 

3.2 The chiral anomaly: external currents 

We now want to study the anomaly of the chiral symmetry when the flavor currents are 
coupled to external sources, i.e. study the WZ term as an action for the flavor brane 
world- volume gauge fields (Al,Ar) which couple in the boundary to combinations of the 
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World-volume field 


T + Tt 


i(T - rt) 


2 


2 


J PC 


0++ 


o-+ 


1" 


1++ 



Table 1: A summary of the world- volume fields, their spins and their transformation properties 
under parity and charge conjugation. In the gauge we will use below (see section 4.1 for more 
details), (Al — Ar) z and (Al+Ar) z are set to zero while the longitudinal component of (Al — Ar)^ 
combines with the pseudoscalars. 



vector and axial currents, (see table [l]). As is usual for Chern-Simons terms, a gauge 
transformation does not leave the action ((T^j) invariant but produces a boundary term. 
As in previous cases, [||, this term is matched with the global anomaly of the dual field 
theory. 

We set the tachyon to its vacuum value ( fUC|) . It is a straightforward although lengthy 



computation to expand the supertrace in (3.6) 
Str e iT 



1 6— form 



^TY e~ T2 { - iF L AF L AF L + iF R AF R AF R + 



+ 2irdr A (A L - A R ) A (f l A F L + X -F L A F R + ~Fr A F L + F R A F^j + 

+ t 2 (A l - A R ) A (A l - A R ) A (F l AF L -F R A Fr) + t 2 (A l - A R ) AF R A {A l - A R ) A F L + 



T 3 dr A (A L - A R ) A {A l - A R ) A {A l - A R ) A (F l + F R ) + 



+ -r 4 A (A L - A R ) A (A L - A R ) A (A l - A R ) A (A l - A R ) A (F L - F R ) 



i 

10' 



-T 5 dr A (A L - A R ) A (A L - A R ) A (A l - A R ) A (A L - A R ) A (A L - A R )} 

(3.12) 

It is now not very hard to find a 5- form ^5 such that dO,^ = Str 1 6 f orm - We can write: 

1 2 f 1 i 

Q 5 = -TV e~ T { -%A L AF L AF L + -A l A A L A A L A F L + — A L A A L A A L A A L A A L + 
o [ 2 10 

1 i 

+ iA R A F R A F R — -Ar AArAArAFr- —Ar A Ar A Ar A Ar A A R + 



iA L A Fr A Fr — iA R AF L AF L + -(A L - Ar) A (F L A Fr + Fr A F L ) + 



1 1 i 

-A L AA l AA l AF l - -A r A A R A A R AF R + —A L A A L A A L A A L A A L + 



- —Ar A Ar A Ar A Ar A Ar 



+ 



+ it 6 dr A 
i 



(A L AAr-ArA A l ) A (F L + F R ) + iA L AA L AA L A A r + 



Al A Ar A Al A Ar + iA L A Ar A Ar A A 



R 



+ 



+ 20 ( Al ~ Ar ^ A ( Al ~ Ar) A ^ Al ~ Ar ^ A ^ Al ~ Ar ^ A ^ Al ~ Ar) 



(3.13) 
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We now consider a gauge transformation of ( 3.13 ). As we argued above, the tachyon 
vacuum configuration ( |2.10| ) we are considering breaks the flavor invariance to the diagonal 
(vectorial) subgroup of U(Nf)i, x U(Nj)r. Therefore, taking Vl = Vr = V, the gauge 
transformation of the fields involved in ( |3. 13 ) reads (see appendix |A] for more details) 



SaA l = 


—i DA 


SaA r = 


—iDA 


5 A F L = 




SaFr = 


[A,F fl ] 



-idA - A L A + AA L 
-i dA - A R A + AA R 



(3.14) 



Inserting (3.14) in ( p. 13 ) we find: 

5 A n 5 = ^Tr e- T ' 2 dA A { (1 + r 2 ) [ - F L A F L - % -{A L A A L A F L + A L A F L A A L + 



+ Fl A Al A A L ) + -A L A A L A A L A A L + F R A F R + % -{A R A A R A Fr+ 

+ -Aii A Fr A j4k + Fr A Ar A Ar) - ^Ar AArA Ar A Ar] + 

+ r 3 dr A [(A L - A R ) A (F L + F R ) + (F L + Fr) A (A L - A R ) + i(A L A A L A A L + 

-AlAAlAAr-ArAAlAAl + AlAArAAr + ArAArAAl-ArAArA Ar)] } 

(3.15) 



As expected, it is straightforward to show that the 5-form 8a&5 is closed 

d(5 A n 5 ) = o 



(3.16) 



The result ( gig ) ensures that the gauge variation of the 5-dimensional CS action (BJ) is 
given by 4-dimensional boundary terms. The UV boundary term (where r — > 0) repro- 
duces the field theory global anomaly 0, as we will see below. However, if we consider a 
background dual to a confining gauge theory, the gauge variation of the CS action ( ft.7| ) 
may receive a contribution also from a potential IR boundary. Because of the confinement 
property, the space-time ends smoothly in the IR, and the brane-antibrane world-volume 
must end before or at the IR end of space. The contribution to the gauge variation of the 
CS action coming from this potential IR boundary would give rise to an additional con- 
tribution to the gauge anomaly, and make the bulk theory inconsistent. It must therefore 
vanish. One is led then to the consistency condition r — > oo at the IR, which leaves us 
with the UV contribution alone: 



SaScs 



iN r 



(2vr) 5 a /3 



24vr 2 



Mink 4 



Tr A(t] L - vr) 



(3.17) 



where we have used the fact that the tachyon vanishes on the UV boundary and reinserted 
a (27TQ;') 3 factor in (|3.15j ). We have defined the 4-form r] as 



i 1 
F A F + - (A A A A F + A A F A A + F A A A A) — —A A A A A A A 



(3.18) 
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Notice that the computation is rather similar to the one in but here we have included 
the dependence on the tachyon . 

In fact, we have checked the anomalous behavior under a transformation within the 
U{Nf)v symmetry preserved by the vacuum. It is possible to find the anomaly in the full 
U{Nf)i x U{Nf)n chiral symmetry before it is spontaneously broken. For that purpose, 
we must consider a more general ansatz in which T is proportional to a unitary matrix 
TT* = T*T = T 2 TLN f - With this milder assumption, the generalization of (3.12) reads: 



Str | 6 _ form = lTre~ TtT { - iF L A F L A F L + iF R A F R A F R + 



+ DT ] ADT AF L AF L + DT ] A F R A DT A F L - DT A DT ] A F R A F R + 



i i (3-19) 

+ -DT ] A A DT ] A DT A F L - -DT A DT* A DT A DT ] A F R + 

- —DT ] A DT A DT ] A DT A DT ] A Dt\ 
60 / 

One should now follow the same procedure as above and generalize the expressions ( |3.13 )- 

( 3.15 ). However, it proves a difficult task to find the five-form O5. In any case, the IR 

term will again vanish when r — ► 00 due to the overall exponential factor. Since in the 

UV t — ► 0, one can find the UV contribution even without knowing the full expression for 

<5a^5- The result is of course ( |3.17| ) but now with independent parameters for the left and 

right gauge transformations. 

'./V f 

SaScs = / Tr [ A L??L - A RVR ] (3.20) 

24vH y Mink 4 

Possibly up to a sign, (which is a matter of conventions), ( |3.20 ) reproduces the QCD 
chiral anomaly, as we now check by using the general procedure sketched in Q and further 
developed in ||. The AdS/CFT conjecture @, g states that S d [A] = W[A] where W is 
the generating functional for current (J A = -rp-) correlators in the boundary theory. We 
can equate the gauge variation of both terms in the equality: JaS'cZ^] is given in ( 3.20) ) 
while S A W[A) = jjx$AA A = jd A xD^K A J A ^ = - f d 4 x A A (D^) A . The anomalous 
divergences of the U(Nf)L,R flavor currents in the presence of sources are, thus: 



a M j?#* = (-) fc *^* -iv („,.„) 



(d,jI r ) =H hL ' R ^*K(\* VL>R) 



Here * stands for the four-dimensional Hodge dual, we have defined and h R to be 
and 1 respectively, and the currents have been decomposed in their abelian and non-abelian 
components following the conventions of appendix [A|. Expressions ( [3. 21 ) reproduce the 
known QCD results (see for example |4Q|| ). 

For completeness, we report here the explicit expressions of the traces appearing in 
equations (3.21). They can be easily derived from the definition (3.18) of 77: 



Tr(ry) = (-) 



p+1 (If AF a +N f F u ^ A F u ^ - ^f ab c A a A A b A F c ^j (3.22) 



7 Notice the different definition of the gauge field with respect to M. A = iA , F = iF where A, F are 
the gauge field and field strength used in M. 
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and 

TY(A a 7?) = (-) p+1 



}_ d abc F b A pc + F a A pU(l) _ ly6c a ^6 f\ j\ c f\ pU(l) _ 2A U(1) A pb A ^c^j + 



_1 C / ab d ecd - / ac e d eM + f bc e d ead ) A b AA c AF d + —f bc h f de ,d ahl A b A A c A A d A A e 
8 V / 16 

(3.23) 

3.3 The U(1)a axial symmetry 



The first equation in ( 3.2 1] ) for the anomalous divergence of the abelian flavor currents in 



the presence of external sources is not complete. In the case of the U(1)a symmetry, an 
additional anomaly appears fl23| . 

To study this U{\)a anomaly, we need to know what bulk field couples to G A G 
on the boundary (here G denotes the field strength of the non-abelian color gauge field). 
We consider then a probe Dq color brane. Since we are studying four-dimensional gauge 
theories, the Dq brane wraps a (q — 3)-cycle JC g s of the internal geometry. Typically, for 
a background dual to a confining color gauge theory with no moduli space nor low-energy 
massive scalars, the Dg-brane embedding is stable only at the boundary, in the far UV 
region. 

The action for this probe brane reads: 

Sprobe = ~T q ! cP +1 y e~V-det(P[ ff ]+G) + T q f C A Tr e iG (3.24) 

where P[g] is the pull-back of the background metric on the world-volume S<j+i of the T>q- 
brane, and as before we have set 2ira' = 1. The coupling to GAG we are interested in comes 
from the expansion of the exponential in the WZ term, and is given by -y- J s Cg-3 A 
Tr(G A G). 

Comparison of this term to the usual Yang-Mills action gives the holographic definition 
for the QCD #-angle (recall that p + q = 12) 



}CD 



^ 2 T l2 - P [ 0)- P (3.25) 



where the integral over /Cg_ p is evaluated in the far UV (the boundary). 

To study effects on the dynamics of flavors related to Tr G A G, we will need, therefore, 
to turn on the RR potential Cg_ p in the background, and look at its interaction with the 
flavor branes. Using Hodge duality, this potential has the same degrees of freedom as C p _i. 
The coupling of this potential to the Dp-Dp flavor branes can then be read from the WZ 



action (|3.l| ) 



S wz = -iT p I A Str e^L. (3.26) 



If we take T to be proportional to a unitary matrix, TT^ = t 2 !^ , as we did at the end of 
subsection 3^, it is easy to write from ( |3.2[) an explicit formula for Str e %:F 



2— form' 

Str e^L. = e- T \i{iF L - iF R - DT^ A DT) (3.27) 
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The coupling in ( 3,26| ) modifies the equation of motion of C p _i, or equivalently the 
Bianchi identity for the Hodge dual Cg- P . Indeed, the terms of the overall action involving 
Cp-i are S = ~2^~^ I \ F p a * F p ~ iT P It, p+1 C p~ 1 a Str ^L-form' from whicn we obtain 

dF W -p = d* F p = 2iK 2 ^T p 5g- p (£ p +i) A Str e lJC ] 2 _ form (3.28) 

where -Fio-p is the gauge-invariant (with respect to world-volume gauge transformations) 
(10— p)-form RR field strength, and 5g- p (£ p +i) is the volume form for the space transverse 
to the flavor branes times a <5-function localized at the position of the flavor branes. 
From ([T2§ we can write: 

Fia- P = dC 9 -p + 2i(-)P +1 K 2 m T p 5 9 -p (E p+1 ) A fti (3.29) 

where f2i is defined by dQ\ = Str eijr | 2 _f orm an d can be explicitly obtained from ( |3.27 ): 

Six = e~ r2 Tr (iA L - iA R + (log T - log T^rdr) (3.30) 



As we show in appendix 0, the only gauge transformations under which Q\ transforms non- 
trivially are the U{1)a axial ones. We thus consider an infinitesimal such transformation 
Al = -A r = iaJL Nf : 



8Ux = 2iN f e- T (da - a d(r 2 )) = d [2iN f e- T a) = duj (3.31) 

2 

where we defined ujq = 2iNje~ T a 

Since -Fio-p is gauge invariant, and as we just showed 5\£li = dujQ, it follows from 
(3.29) that the RR-potential cannot be invariant under flavor brane world-volume gauge 
transformations [38]: 

^aCq^p = —2iK 2 w ^T p 5g- p (S p+ i) ljq , (3.32) 

which, inserted in ( 3.25| ) and using the identity 8ir 2 Ti2- p T p k 2 = 1 (in units of 27ra' = 1), 
results in the following transformation relation for the QCD theta angle under U(1)a 
transformations : 

SaOqcd = -i uq\ uv = 2N f a . (3.33) 

Because of the boundary coupling j d A x e ®^2 Tr(G A G), the formula ( |3.21| ) for the di- 
vergence of the U(1)a current J^ 1 '^ = J^ 1 ^ — Jr has to be corrected in order to 
take into account the non-trivial transformation of 9qcd (the U(\)y current Jy M is 
unaffected) : 

Va^" = ^* Tr( - r lL + r 1 R) + ^ e^Tr su(Nc) {G^ 2 G^ A ) (3.34) 



8 The following formula has been computed for the case where all the quarks have the same (possibly 
vanishing) mass, which is the situation considered in this paper. Since in section i.4 we will relate Qi to 
the would-be Goldstone boson r[ , we should expect that in a more general case Tr(logT- logT 1 ") will be 
substituted by an expression of the form log(det T(T t ) _1 ). We will not pursue this computation here. 
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3.4 A holographic view of the Coleman- Witten theorem 



In |22j, Coleman and Witten proved that, under a few assumptions, in N c — > oo massless 
QCD, chiral symmetry is spontaneously broken U(Nf)L x U(Nf)n — > U(Nf)v- 

One of the main results of the previous section is that, if the theory is confining, 
the tachyon has to get a (diverging) vev near the IR even if T — > in the UV. Since 
T transforms in the bifundamental representation of the flavor group, (T) ^ means 



that the chiral symmetry is broken. We showed in section 2.2 that (T) is proportional 
to the identity matrix, and consequently the flavor symmetry is broken down to U{Nf)y. 
Confinement, therefore, implies spontaneous chiral symmetry breaking, reproducing in a 



holographic language the result of [22 



To make the analogy between these two different approaches more evident, we can 
compare the assumptions of p2| with those we made in the present holographic setup. 



In order to do this, we start by reviewing the five assumptions of [22]: 1) The large iV c 
limit of QCD exists. 2) In that limit the theory is confining. 3) There is a single order 
parameter for the breaking that is a quark bilinear and transforms in the bifundamental of 
U{Nf)i x U(Nf)R. 4) Its vev can be found by minimizing some potential. 5) The potential 
is assumed not to have degenerate minima. 

For the holographic setup, assumptions 1) and 2) are also necessary, and of course, 
one must further assume the (non-trivial) fact of the existence of a holographic string dual 
of large iV c QCD which fits in the general framework we have introduced. It is appealing 
that assumption 3) is automatic in the setup of our work since it is a well established fact 
that the lowest state of a brane-antibrane system is the open string tachyon. Concerning 
assumption 4), the string theory dual provides, in principle, a way to determine the vev 
of the quark bilinear: the bulk vev of the tachyon field is determined by solving the bulk 
field equations with the UV boundary condition corresponding to m q = fl2.9|) and the IR 
consistency condition r|//j — > oo. Then, one can read the value of the quark bilinear vev a 
from the UV behavior of t(z) (|2.9| ). This computation requires the DBI part of the action, 
and is similar to that in |3[ It was discussed in section 2.2 . 



In principle, one expects that the UV and IR conditions determine uniquely the bulk 
vev of the tachyon and therefore the (qq) condensate of the dual theory. Nevertheless, 
in the very general setup we are considering, we cannot prove that the solution for the 
(qq) is indeed unique, so this must be taken as a further assumption, analogous to the 
non-degeneracy 5) listed above. 



We remind the reader that in section 2.2 we also made the (reasonable) assumption 



that the stable minimum given by the vacuum (|2lC| ) is the real vacuum of the theory. 

The demonstration presented in this paper, namely that confinement is a sufficient 
condition for spontaneous chiral symmetry breaking, can be thought of as a reformulation of 
the geometrical picture of || . There it was argued that if one places branes and antibranes 
in a background that smoothly ends at some point of the radial coordinate (and therefore 
is confining), they must necessarily recombine in the IR, sparking the breaking of chiral 
symmetry. The advantage of the present formulation stands in the fact that the tachyon, 
the scalar responsible for the symmetry breaking, is explicitly taken into account which, 
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for instance, allows to introduce a bare quark mass. Notice that, if one is in the conformal 
window of N = 1 or J\f = QCD, then there is no IR boundary and it is consistent to 
take a vanishing vev for the tachyon as in [24, 19]. However, as in p^] or |J, our argument 
does not rule out having spontaneous chiral symmetry breaking without confinement. 

3.5 The effects of a non-trivial quark mass 

In [39 1, it was explained how an expression like ( |3.21| ) can be related to the anomalous 
three-point function in the case of N = 4 SYM where the global symmetry associated to 
the current J is SU(A)r. Analogously, in our case one can relate ( |3.21| ) to the anomalous 
massless QCD three-point function. This anomaly equation in the large N c limit was used 



by Coleman and Witten [22 1 to prove the existence of massless Goldstone bosons, and 
therefore the spontaneous breakdown of the symmetry U(Nf)L x U(Nf)n — > U(Nf)y- 
Clearly, if there is a non-zero bare mass for the quarks, there cannot be spontaneous 
breaking nor massless Goldstone bosons, so (3.21) must be modified when m q ^ 0. In 



section |4.3| , we will show by analyzing the equations that define the mesonic spectrum that 
there are massless Goldstones if and only if m q = 0. The goal of this section is to describe 



holographically the modification of ( [3.21 ). 



The axial current is classically non-conserved when m q ^ 0. In the holographic picture 
this shows up because a gauge transformation of the bulk gauge fields is accompanied by 
a gauge transformation of the tachyon ( |A.6| ), which contributes to ( p. 21 ) iff m q ^ 0. In 



the rest of this section we make these statements more precise, assuming that the space is 
asymptotically AdS: 

9xx (z) ~ R 2 Ad s/ z 2 + ... g zz (z) ^ R 2 AdS l z 2 + ... (3.35) 

The antihermitian part of the tachyon matrix is dual to (— i times) the pseu- 

doscalar current, see table 1. We consider a perturbation around the vacuum T = re %d 
where 9 is proportional to the (hermitian) pion matrix, 9 oc J^ N II + 7r a A a . At first order, 



T ~2 = ir9. There must exist a boundary coupling contributing in the W generating func- 
tional J d 4 x (—i<fi A J A ) z=0 , where J5 is the pseudoscalar (pion) current J5 = iqj 5 q+iq r y 5 X a q 
and 0o is proportional to ir9. As argued in section we require that the mass of the 
tachyon, at least near the boundary, is Tn^,R AdS = —3 since qq and iq^q are operators 

of dimension 3. Following Witten's prescription [||, one has to couple to the pseudoscalar 

3 

current in the boundary, 4>q = R\dSz^ Tcan ^i since f° r a scalar of the above cited mass, the 
non-normalizable behavior is ~ z. Inserting the value of r can given in (|2.9|), the contribution 
to the generating functional is: 

'' d 4 x (m q 9 A j£) (3.36) 



We now compute how this expression transforms under an axial gauge transformation Al = 
—Ar = A = iA A X A = iaR + iA a \ a . From ( |A.6| ), the infinitesimal gauge transformation 



V = e of the tachyon around the vacuum ( 2.10 ) is given by d\T = ArT — TAi 



t(Ar-Al) = — 2At, which in turn can be expressed as a transformation of the pion matrix 9 
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as: 5\0 A = (2iA) A = -2A A , which, inserted in (|3~36|) , yields / d L 
This leads to a modification of (3.21) and ( 3.34| ): 



x m„ 



2aJ 5 f/(1) -2A a J 5 a ). 



df,J A 



A 

D^aT 



U(l) 



-2m q J 5 



-2m q J% + 



+ 



24vr 2 



* Tr (jil + rj R ) + 



16vr 2 



eMl " ' M4 Tr 5(7 ( N c ) ( G> x ^ 2 G w w ) 



24vr 2 



*TV(A a (r ?L + r 7R )) 



(3.37) 



One can now relate these expressions to the three-point function in the presence of a non- 
trivial quark mass. The term with m q in ( |3.37| ) results in an extra term for the three-point 
function. The presence of this new term invalidates the argument that Coleman and Witten 
used in the m q = case to show that the three-point function has a pole at zero momentum. 
Thus, as expected, the three-point function can be analytic at zero momentum, consistent 
with the absence of massless Goldstone bosons. 



4. Features of the mesonic mass spectrum 
4.1 Equations determining the spectrum 

We now study small fluctuations of the different fields around the vacuum configuration. 
They correspond to the mesons of the field theory. We will just consider quadratic terms in 
the action, which are enough to find the mass spectrum, but not the couplings. At this level, 
the non-abelian action ( |2.4| ) is just the sum of Nj copies of the abelian one (|2.l| ). Thus, 
we use (2.1) in the following and there are Nj copies of each of the towers of resonances 
that we will study. 

We start by considering excitations for the fields 6,A L ,A R around the vacuum of 
section 2.2. Fluctuations of the modulus of the tachyon r will be briefly mentioned below. 
We expand the simplified version of ( |2.1| ) (neglecting the scalars Y/^ and assuming there 
is no -B-field). We define: 



9zz + -(d z r) 

TV 



e~^ = e~^ffV(r 2 ), g zz 
In order to use a notation similar to the one in |T3|, we choose a gauge: 



(4.1) 



A 



R 







(4.2) 



and define: 



V M 



aL i aR 
M + A M 



A 



M 



aL _ aR 

M M 



2V2 



(4.3) 



and V^ V ,A^ U as the (abelian) field strengths of V^A^. Expanding (|2.lD , we find: 



S 



d xdze 



~9UV^ U V^ + + g xx gj ({dzVrf + (dzA^) 2 ) + 



\g xx ~giv 2 {d„e + 2A^f + \gl x ~gjv\d z ef 



(4.4) 
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Vector sector 

The vector (1 ) sector is decoupled from the rest in the action (|4.4j ). The equation of 
motion for can be solved by expanding in modes: 



^ = E^) v i n V) 



(4.5) 



with: 



d z (e ^9xxgzz 2 d z i> {n) ) +m 2 n e *g2 z ip( n ) = (4.6) 
The physical modes are those which yield a finite action when integrating in z. Substituting 



(4.7) 



|4.5| ) in (4.4), one finds a tower of massive vectors: 



provided the normalization conditions: 



. 9 



ZIR 1 - l 

dz-e-*ghtf n) 



ZIR 



dze ^gxxgzz'idztpin)) 2 = m 2 n 



(4- 



Integrating by parts the second expression and using (|4.6|) , we consistently obtain the first 
line of flO]). 

Axial vector sector 

The vector field fluctuation can be split in a transverse and a longitudinal part, = 
Ah + Af,,, with d^Aj^ = 0. We first consider the transverse part, corresponding to 1 ++ 
excitations. We expand it as: 



4 = E A w V) 

n 

which results in a tower of massive axial vectors: 

s = - /^E ll^K) + <^ {n) 

subject to the normalization conditions: 



(4.9) 



(4.10) 



ZIR 



dz 



e-^gxxgzzHd z Af n) ) 2 + e^g xx glv 2 Af n) 2 



1 

2 



(4.11) 



We have chosen this normalization such that when the fields are promoted to their non-abelian gen- 



eralization Vjr — > Vp"' a \ a , the components V^ n; a are normalized in the standard way, since we are using 
TrA a A 6 = \8 ab . We apply this prescription to all the different towers of modes. 



;(«) a 
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and the second order differential equation: 



d z (e-*g xx gJd z Af n) ) + mle~*ghAf n) - e^g xx glv 2 Af n) 







(4.12) 



It is interesting to notice at this point a difference with respect to constructions like the 
Sakai-Sugimoto model J7|. In that kind of models, the vector and axial vector mesons are 
solutions of a single second order differential equation with different matching conditions 
at the IR. In this case, they satisfy different equations 



( 4.12 ). If the chiral sym- 
metry were not broken (v = 0), the equations would become degenerate and vector and 
axial vector mesons would have the same mass spectrum due to the unbroken symmetry. 
The differential equations found in this construction coincide with those of the AdS/QCD 
models (l^, [l4|] up to the z-dependent dilaton, metric factors and tachyon vev, which we 
have kept general. 

Pseudoscalar sector 

The modes coming from 9 and the longitudinal part of combine to give a single tower of 
resonances with the quantum numbers of pions l ". The equations of motion from ( |4.4| ) 
can be solved by expanding: 

4 = -E^n)W(« (n) K)) 
n 

9 = 2]Ttf (n) (z)a("V) (4.13) 
n 

where the functions (fr n \ and i9( n ) satisfy the following coupled differential equations 



- i ~ i 

d z (e' 4 'g xx g Z z ! <9^( n )) + ^ g xx gi z v 2 {"9^ n) - tp( n) ) = 

gxxV 2 d z $ <yn) - m 2 n d z (p( n) = 



Inserting (4.13) in 

l * z IR 



dz 



o 



4|) and normalizing: 

- _i - i 

e'^gxxQzz 1 {d z f( n )) 2 + e'^gxxgizv 2 ^^) - V3 (n )) 2 



dz 







*glx~gJv 2 {d z d (ri) f 



one finds a tower of pseudoscalars in the four-dimensional theory: 

S - - I d 4 xY, [(^("l) 2 |m^(a("l) 2 



(4.14) 

(4.15) 
(4.16) 

(4.17) 



Scalar sector 

The scalar mesons (0 ++ ) come from excitations around the vacuum of the form 5t 
S(x^,z). Expanding ( |2.6| ) up to quadratic order in S, one finds an action: 



S = -2 I d 4 xdze 



ld 2 V 



2 V 



+ -gl x gzz 2 {gz 

7T 



T=v/2 9x*9zzb + 2 ^ y 
~ 1 1 



g 2 xx fJ{d z v)S(d z S) + 



t=v/2 



-{d z vY){d z sy + -g xx ~giz{d„S) 



(4.18) 
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from which one can straightforwardly extract the linear equation of motion and normal- 
izability condition. Since we will not need them in the following, we do not present them 
explicitly. 

4.2 On "linear confinement" (for highly excited mesons) 

In a theory with linear confinement such as QCD, one expects that the squared masses 
of highly excited hadrons grow linearly with the excitation number n: cx n. This 



behavior is difficult to find from holographic models, that typically yield tx n 2 [41]. 
In jl(| it was shown that, by appropriately tuning the behavior of the metric and dilaton 
in the IR, one can indeed find the expected behavior. Moreover, it was conjectured that 
such behavior could come from closed string tachyon condensation, but no concrete model 
producing such IR asymptotics of metric and dilaton has been found. 

We will show here that the relation m 2 cx n is automatic in our construction. Therefore, 
the linear confinement relation for mesons comes from open string tachyon condensation 
on the flavor brane world-volume. In this section we just state the ideas and results while 
details are relegated to appendix [E|. 

The argument relies on the fact that the meson excitations feel the effective open string 
dilaton and metric (see eq. (|4.1| )) rather than just the closed string dilaton and metric. In 
the IR, t — > oo so g zz = g zz + ^(8 z t) 2 f (3 2 r) 2 even if g zz has a single pole. Moreover, 
9xx goes to a constant, which, on general grounds, can be identified with the QCD string 
tension. Reinserting a' and defining a new radial variable u = \j ■ k Tqcd T ^ z ^'' we nave ' near 
the IR, where u, r — > oo: 

ds 2 = 2ira' (r Q cDdxl 3 + ^(d z T) 2 dzA = 2-koITqcd {dx\ z + du 2 ) 

e~* ~ V(r 2 ) ~ e^ 2 ~ e "^« 2 ^ ^R u ^ (4.19) 

It turns out that the natural radial variable in the IR is proportional to r. The quadratic 
growth of the dilaton ( |4.19| ) is the IR behavior advocated in [16] to account for linear 



confinement. In fact, using the WKB approximation to compute the masses of highly 
excited bound states one obtains, for the vector meson tower (see appendix [E]): 

m 2 w 2-kTqcd n (large n, vector mesons) (4.20) 

This is a quite general result that only relies on having a confining background and 
therefore diverging tachyon, on the large r gaussian behavior of the tachyon potential and 
on the DBI action. We find very encouraging that our construction naturally accounts for 
the physically expected relation of linear confinement, including the correct multiplicative 
factor. 

However, there is an important caveat: if one repeats the same computation for axial 
vector mesons, eq. (4.12), one finds: 



16, 



m n = 2iryl H — ^Tqcd n (large n, axial vectors) (4-21) 
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This cannot be physically correct since vector and axial mesons should asymptotically 
have equal masses at large n, as chiral symmetry is restored for excited hadrons G^]. The 
problem may arise from the fact that we are using the DBI action outside its range of 
validity. This point definitely deserves a better understanding. 

Finally, it is important to stress that in view of the confining IR behavior, there is 
no need to impose an arbitrary IR boundary condition. The IR condition for physical 
excitations is simply the normalizability of the action [16]. 

4.3 Goldstone bosons 

Large N c QCD with massless quarks has a set of Nj massless pseudoscalars which are the 
Goldstone bosons of the spontaneous breaking of the U{Nj)a- We will generically call 
them pions. In this section we show how they appear in this construction 10 . They are the 
solutions to equations ( [4.14| )-( |4Tl~6| ) with m n = 0. Notice that, as expected, there cannot 
be massless (axial) vectors due to (p~g|), ( |4. 1 1| ) . 

We generalize the analysis in ( [4.13| )-( [4.17| ) with m n = to the non-abelian Nf > 1 
case. Define the pion matrix as the generalization of the a^°\x u ) in ( [4.13| ): 

7T (xn = ^=^ + 7T a (xn\ a (4.22) 

To lowest order, the action for the pions is just the one for a set of massless scalars: 

S = - J d 4 xTi (d^Tr) = - J d 4 x Q^T/dV + ^d^d^nA (4.23) 



Due to Q4.14 ), i?( ) must be a constant. We will shortly see that it is related to the pion 



decay constant. We define a quantity: 

£(z) =¥>«)) (*) (4-24) 
which, regarding (|4.14| ), ( 4.15|) , satisfies: 



l 9xx~g z hzil{z)) - e-*g xx giv 2 az) = , (4.25) 
-e-^g xx ~ 9 ;h(z)d z ^)\ z=0 = 1 (4-26) 



For the second line, we have integrated by parts (4.15) and taken into account that for £(z) 
we must choose the normalizable mode in the IR. 

The question of whether there are massless modes boils down to the existence of a 
function £(z) satisfying (|42g),(B). From flO]), it is obvious that the answer depends 



on the UV behavior of £(z), which, due to ( |4,25| ) depends on the UV behavior of 

We now give a heuristic argument suggesting that if quarks are massless, there are 
indeed massless Goldstone bosons. Since r is dual to qq, the UV behavior of r is related to 
the quark mass and condensate. If there is no quark mass, r is normalizable, i. e. it vanishes 



One of the iV* pseudoscalars is massless only in the strict N c — > oo limit. This is described in section 4.4 
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fast enough in the UV. It is natural to think then that the second term in ( 4.25 ) can be 
neglected near the UV, allowing for a solution to ( |4.25| ) and ( |4.26| ), which is, asymptotically, 



£(z) = -0(o)+0p)/(*) 



(4.27) 



The constant i?( ) has to be determined by the IR condition, we have defined f{z) such 

_ i 

that d z f(z) = g X x9zz )~ X and used the residual gauge invariance to fix y>(o)(0) = 
and therefore /(0) = 0. If the quarks are, instead, massive, v is larger in the UV and the 
second term in ( [4.25 ) cannot be neglected, which hinders the existence of a solution as the 
one just described corresponding to a massless Goldstone. 

Just as a clarifying example, we consider a five-dimensional space which is asymptot- 
ically AdS. In this case, all expressions asymptote in the UV to those of |l3|, where one 
can explicitly check that the heuristic reasoning above is valid. 

The pion decay constant 

Consider m q = so that there are indeed massless pions. The pion decay constant can be 
related to the pole of the axial current two-point function at zero momentum: 



nA(? 2 ) = £ 



f 2 



J n 

q 2 +M'i n ' q 2 



+ 



(4.28) 



One can compute this correlator following the AdS/CFT prescription, as in |l4}| . This 
is done by computing the on-shell action giving the appropriate boundary condition to the 
field dual to the axial current, i.e A^. When q 2 = 0, the bulk equation for Aj^(z) is ( 4.12j ): 



d z {e-*g xx ~gJd z Af n) ) - g xx ~gl z v 2 A\ n) = (4.29) 
Substituting in (4.4) and deriving twice with respect to (as defined in ( (4.9| )), one gets: 

fl = lim (-e-*g xx g7M ± (z)d z A ± (z)\ z=e } (4.30) 
where A (z) is a solution to eq ( [4.29| ) subject to the UV condition A | z=e = 1 and the IR 



normalizability condition. Notice there is an extra factor of 1/2 in (4.30) with respect to 
what one would get from (|4.4|) . This is again because we use conventions suitable for the 
non-abelian generalizations of the fields, see footnote |9[ 

Two remarks are in order: even if the f n in ( 4.30j ) depends on quantities at z = 0, the 
value of Sn depends on the full metric, since one has to select the well behaved IR mode. 
This typically involves a non-trivial numerical integration. Notice also that f 2 is of order 
N c since inside there is a D-brane tension which scales as gj 1 ~ N c . 

Comparing Q4.25|) , (|4.26|) , (|4.29|) , flOol ) and using the fact that both £(z) and A ± (z) 
must follow the normalizable behavior in the IR, it is straightforward to conclude that up 
to an unimportant sign, £(z) = — f^ 1 A- L {z). This equality at z = yields: 



(4.31) 
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The Gell-Mann-Oakes-Renner relation 



We now show how, considering asymptotically AdS, one can obtain the GOR relation [43], 
which gives the masses of the pions when the quark mass is small but non- vanishing. The 
argument is very similar to the one in [13] but details are different. In the following we 



assume m q ^ but only keep terms linear in m q . 

We have to solve (4.14)-(4.16) for small m n = m^. This solution can be obtained as a 



perturbation of the m n = case studied at the beginning of this section. As opposed to 
that case, "d cannot be a constant and, from ( 4.16| ), one can show that ?9| z= o = is needed. 



Then, using fl4.14|) 



*{z) = ml I dz^4 = 4^ f dz- ^= (4.32) 



o 



9xxV 2 2 J (m q z + az 3 ) 2 

For the second equation we have used that the integral is dominated by the small z region 
so we can substitute the asymptotic value of the different functions. To obtain d z ip(z) we 
have used its value in ( |4.27| ) and for v 2 we have substituted ( |4.3[ ) and ( |2.9| ) . As a consistency 
check, notice that the region where the integral above has significant support is around 
z ~ y^niq/a, so taking m q small enough, the integral only probes the asymptotically AdS 
region. 



For z 3> y/m q /a, the function $(z) goes to a constant which has to be the one of the 
massless case ( |4 . 3 1| ) . Putting everything together, we get the known expression: 



9 4m„(j 2m a (qq) 

J it Jn 



where we have substituted ( B.10| ). 
4.4 The mass of the rf 

We now return to the m q = case, where, in principle there are N 2 massless Goldstones. 
In section |3.3| it was shown how our setup correctly reproduces the O (jf^ anomaly of 
the U(1)a axial symmetry This anomaly implies that the (generalization of the) rj' 
meson, the would-be Goldstone boson corresponding to the diagonal U{1)a subgroup of 
the spontaneously broken U(Nf)A, has a mass of order j^-, and is therefore massless only 
in the strict N c — > oo limit. 

The r/ mass appears in the present holographic setup via a Stuckelberg mechanism. 
The reasoning we will follow is very similar to the one of Q (see H, ^4| for related work 
and [^] for a recent discussion in deconfined theories with broken chiral symmetry). 

The Cg-p can only appear in the action in the gauge invariant combination ( 3.29| ), 



Integrating this expression in the space composed of the 9 — p non-Minkowski dimensions 
that the color branes wrap plus the radial direction we obtain: 

(27r) 2 T 12 _ p [ Fio_ p = Oqcd + i(-) p+1 [ (ni),dz = 9 QCD + R p ^%' • (4.34) 



We have used ( 3.25|) and J M dCg- p = J^- Cg_ p where the last integral is evaluated 



in the UV. For a 3+1 confining theory defined on Dg-branes with q > 4 as in |33L Eq] this 
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equality holds automatically because the space closes off smoothly at the IR and /Cg_ p 
is the boundary of A^io-p- On the other hand, if one wants to build some kind of five- 
dimensional model in the spirit of [13, 14, 15, 16, 17|, one would have to impose that the 
RR-potential vanishes in the IR boundary to find the same condition. 

The integral in ( [4.34| ) is computed substituting T = re = re i2, M :r '') — Te l u in 



(|3T30[) : 

(Q 1 ) g dz = J e^Tr ((logT- log T^rdrj = -^-Tr(vr) J e^rdr 



— V 



(4.35) 

The contribution of Oqcd to the vacuum energy density can be computed by integrating 
the kinetic term of the RR Cg-p-form [47], and must appear through the gauge invariant 
combination ( [4.34 ). Thus: 



2 



S = -\ J d 4 x(9 Q cD + (-T^^v'] (4-36) 



where the topological susceptibility of the vacuum is, by definition, the second derivative 

dffj 



of the vacuum energy density in the glue theory without flavors x = ( j£ £ ) ■ ^ 



no quarks 



QCD _ 

can be computed from the supergravity action []47|]. The expression ( 4.36| ) reproduces the 
Veneziano-Witten formula for the mass of the rf [48] 11 : 



, 2JV f 

m v' = -fTX (4-37) 

J TV 

5. Comments and discussion 

Compared to the models Q (Q) where flavors where introduced with D7 (D6)-branes in 
the background of D3 (D4)-branes, our construction has the advantage that the breaking 
of the full non-abelian chiral symmetry is described. In fact, one may think of the models 
presented in those papers as the result of a Dp -Dp system in which the tachyon has fully 
condensed leaving a D(p — 2) brane as a vortex. Since a non-trivial tachyon breaks the 
U(Nf)A, this is not present in the final configuration and only a remnant U(1)a is left as 
the rotation symmetry around the vortex (for a related discussion, see ]|l8|]). 

The model of Sakai and Sugimoto |7], ||] and its generalizations share many properties 
with our setup, but there are also important differences. In [0, §, it does not seem possible 
to include a bare mass for the quark, so clearly there are aspects of the QCD symmetry 
breaking that are not present in that setup. Another phenomenon that we have shown to 
be intimately related to the physics of the tachyon is the existence of massive resonances 
with the quantum numbers of pions '". These are not present in the model of 0, ||] . A 
heuristic way of understanding the differences is that the description of [0, H is done already 
in the broken phase (with the brane-antibrane pair reconnected), so that information about 



11 The difference in a factor of 2 with respect to is due to a factor \/2 in the definition of f n . In our 
conventions (0| J ^' = ^/2A f /p M /,,/ and f n i = up to suppressed O^N' 1 ) corrections. 
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the order parameter is lost. It would be interesting to try to include the tachyon in such 
type of models, generalizing the formalism of this paper to non-overlapping brane-antibrane 
with excited transverse scalars. For some progress in this direction, see [p7Tf| . 

In a sense, the model presented in this paper can be thought of as a string theory 



construction in which the AdS/QCD models 14 1 are embedded. We have seen that the 
five-dimensional field spectrum is the same. Expanding the square roots for small values 
of Am and T, one gets the same kind of action as in [O, 14]. In fact, the equations to 



determine the spectrum (|4.6f ) , ( |4,12|) , (|4.14|) reduce to those in [[D|, |J] if one considers AdS 
metric and constant dilaton. On the other hand, the equation for the scalar mesons ( 4.1$) 
is different from its analog in |49j] . Small T, as argued above, corresponds to the UV. The 
successful features of |l3|, |l4|] come from the UV. We therefore conclude that a model built 
along the lines described in this paper in asymptotically AdS space can capture the good 
physical features of [13, 14] while the IR arbitrariness of those models is lifted due to the 
condition t\ir — > oo. This condition fixes the quark condensate in terms of the quark mass 
and removes the extra non-physical parameter appearing in 13, 14]. Moreover, the brane 
physics automatically provides the 5d Chern-Simons term. 

By applying the formalism described in this paper to a concrete model (i.e. to par- 
ticular expressions of the metric and dilaton), one could find numerical estimates of QCD 
observables such as meson masses or couplings. Nevertheless, as stressed several times, 
enforcing a well behaved asymptotically AdS space, implies rnx^AdS = — ^ because the 
quark bilinear has dimension 3. Since = —-^j, we need R 2 AdS = 6c/ and the space has 
large curvature, of the order of the string scale. The complete dual background, then, is 
not a solution to just Einstein equations, but higher derivative corrections to the equations 
of motion have to be introduced. This is a very difficult, if at all possible, task. 

The drawback of large curvature implies that meson masses or couplings which would 
be numerical results obtained from the DBI action cannot be considered trustworthy. This 
same problem has arisen when trying to build non-critical holographic models. Neverthe- 
less, Einstein-like equations have been used to extract qualitatively correct results |24|, |T0| , 
5C|. Considering the impressive quantitative precision of predictions in models like |8|, 



13| , 14 1, it would be interesting to build at least phenomenological models incorporating 



the tachyon physics. Constructions like [17|] could be a starting point. 



In any case, apart from quantitative computations, we expect the general features 
derived from the DBI action to hold. Moreover, since the WZ term is topological, the 
results derived from its analysis in section |3| hold even if the curvature is large. 

We end by commenting on some additional open problems. First of all, it would be of 
major interest to add finite temperature and describe the deconfinement phase transition. 
The physics of fundamental matter in this regime is very rich and has been studied holo- 
graphically (see [ pl| ] and references therein). It may be possible to generalize the analysis 
to the kind of setup described in this paper. It would also be nice to understand the physics 
of rotating strings and Wilson loops. Another possible generalization is to introduce a large 
number of flavors Nj ~ A^ c , to go beyond the quenched approximation. This may be done 
along the lines of [24, 19, 52]. In fact, the general form of our expressions should guarantee 
that they can account for any such backreaction, as long as the theory with fundamentals 
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is still confining. Finally, it would also be interesting to make contact with holographic 
M = 1 theories built on the cigar, where some exact world-sheet computations can be done 
despite the large curvature, pS]] . 
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APPENDIX 
A. Conventions 

For the U(Nf) generators X A , A = 0, . . . Nj — 1, we take 



(A A )t = \ A Tr(X A X B ) = X -b AB (A.l) 



This in particular fixes the normalization of the J7(l) generator: A = — ^==1, where II is 

the Nf x Nj identity matrix. For the non-abelian SU(Nf) generators A a , a = 1, . . . , Nj — 1, 
we have 

[A a , X b ] = if ab c X c Tr(A a {A b , A c }) = d abc (A.2) 

where f ab c and d abc are, respectively, the structure constants and the normalized anomaly 
Casimir for SU(Nf). Because of ( |A.1| ) f ab c and d abc are real numbers. 
We define gauge fields to be hermitian 

= A A X A = A^JL + A^X a (A.3) 

In differential form notation, the field strength and covariant derivative read then 

F = dA — iA A A D = d-iA- (AA) 

where A- indicates the representation-dependent action of the gauge algebra. In particular 
the Bianchi identity reads DF = 0, and the covariant derivative of the tachyon is given by 

DT = dT + i TA L - iA R T DT ] = dT ] - iA L T ] + i T ] A R (A.5) 

Under gauge transformations, the left and right gauge potentials, and tachyon trans- 
form in the following way 

A L ^V L A L vl-idV L vl , A R ^V R A R vl-idV R V^ 

(A.6) 

t^v r tvI , rt-ybrtyj , v L vl = v R v^ = i 

An infinitesimal gauge transformation is defined as V e (x) = e eA ^ ~ 1 + eA(x) and the 
gauge transformation of a field as A — ► A + cSaA. From (|A.6| ) we have then: 

5 A A = —i DA = -i dA + [A, A] 
S\F = [A, F] 

1 J (A.7) 
5 Al T = -TA L 

5 Ar T = A R T 

Notice that the generators of gauge transformations are antihermitian. When we decom- 
pose them in their U(l) and SU(Nf) parts, we will write then 

A = iA A X A = ial + iA a X a (A.8) 
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with a and A a real parameters. In particular we have, from (|A.3|) and (|A.8| ) 

5A U ^ = d^a and SA% = {D fl A) a (A.9) 
Currents are defined to be hermitian. We decompose a U(Nf) flavor current as 

J M = J°A° + J*A a . (A.10) 
This decomposition corresponds to defining the A t]l component as 

Jl,Rh = ^coiors(iq^^^X A q). (A.ll) 
For the U(l) component J°, this would be a strange normalization 

JU-^lW-ii^). (A.12) 
We therefore define a rescaled J^ 1 ) = y^2Nf,J°, such that J now reads 

(A.13) 

Notice that the normalization of A u ^ in ( A. 3 ) has been chosen in such a way that the 



j — 1 jU^if + j a \ a 
M 2iV/ M 



boundary coupling of the current to the gauge field reads 

2 y c/ 4 xTr(J^A M ) = J d 4 x (jVMpaVW + J a »A a ^j (A.14) 

B. Determination of (qq) from holographic renormalization 

When solving the equation for the modulus of the tachyon in asymptotically AdS space, 
we found that it depends on two integration constants m q and a, see Q2.9D . The constant 
associated to the non-normalizable mode m q can be immediately identified with the quark 
mass. On the other hand, a is related to the quark condensate (qq) in, in principle, a non- 
trivial way. Schematically, (qq) = — where S denotes the on-shell action. However, 
the on-shell action is UV divergent, and must be renormalized by adding covariant coun- 
terterms. This is done by following the so-called holographic renormalization procedure 



(for a review, see [|53|). In the following, we adapt to our case the method of [54], where 
holographic renormalization was applied to probe flavor branes. In fact, it will be enough 
for our purposes to consider a simple case in which the scalar r does not depend on the 
Minkowski x^-coordinates. Also for simplicity, we consider that the metric and dilaton only 
depart from their AdS values at an order that does not contribute to UV divergences, i.e. 
9xx (z) = R 2 AdS /zHl + 0(z 5 )), g zz (z) = R 2 AdS /z 2 (l + 0(z 5 )) and e~^f = e -*>(l + 0(z 5 )). 

We write the action in terms of a canonically normalized, rescaled tachyon, as defined 
in (P): 

T = C T can , (? = -^T (B.l) 
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The action (with a UV cutoff e) in terms of these quantities reads: 



S, 



ttR 



AdS 



reg 



zir 



1 



d x I dz—^e 



1 + -c 2 z 2 (d z T ( 



2c 2 ./ 7, ~"~ z 5 ~ V ' 3' 

In the asymptotically AdS region, the equation of motion for r can reads: 



(B.2) 



-3r rj 



z 5 d z (z 3 d z T can ) + — z 



(d z T can )d z (z 2 (d z T can ) 2 ) 



3 l+ 2 C 2 Z 2 (d z T can ) 2 



and can be solved in series for small z as 

1 



Tr, 



z 



$ (0) + o C Z lQ g Z$ (0) + Z $ (2) + 



(B.3) 



(B.4) 



Inserting back this result in the on-shell action, one finds UV divergencies as e — > 0. They 
have to be subtracted by adding the following counterterms localized at a z = e slice: 

a nRAdS f ,4 / 



Scti 

5*ct2 



8c 2 

ItR-AdS f ,4 /^— 2 

CI Xy 7 ''"can 



6 



nRAdSC 2 f ,4 



18 



y (i 4 x^7 (log e) r, 



4 

can 



(B.5) 



where 7 is the determinant of the induced metric in the z = e slice, i.e. \/— 7 = R\dS e 4- 

(B.6) 



As pointed out in |54|| , one can also add a finite counterterm: 

Set?, = I d 4 xa ^7 r t 



r 4 

' can 



where a is some constant. It was argued in []54| that different values of a correspond to 
different renormalization schemes. Defining: 



S su b — S re „ + SctO + S c tl + S c t2 + s, 



ct3 



(B.7) 



the quark condensate is given by: 
(qq) = lim 



AdS 5r can (e) 



■|^L*« + 4 d ^?0)(^^ - 4«) (B- 



Noticing from Q2.9|) , ( |B.4| ) that $( ) = ^-Ads m< i'i ^(2) = ^Ads a an< ^ substituting -ff^ d5 
6a' = r, we find: 



7T 



(gg) = -2a + ^m 3 q (c 2 x / '- - 4a) 



(B. 



In section 4.3 we have used the value of the condensate for small m q which is unambiguously 
given by: 

(qq)^-2a, (m, - 0) (B.10) 
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C. Comments on the superconnection formalism 



In section || we have worked with the WZ world- volume action for stacks of brane-antibrane 
pairs and used a superconnection formalism which makes the notation quite compact. 
In this appendix we review some definitions and properties of this construction. The 

supermatrices are 2Nf x 2Nf matrices of differential forms M = ( ^ ^ \ where the blocks 

A, B, C, D are Nf x Nf matrices. We deal with two kinds of supermatrices: if the blocks 
in the diagonal A, D are composed of odd (even) differential forms, then, the off-diagonal 
ones B, C consist of even (odd) forms. For instance, A (J-) defined in ( |3.2[) are matrices 
of each kind. 



The multiplication of supermatrices is defined as [36]: 




AA' + {-f'BC AB' + {-) d 'BD' ' 
DC' + (-) a 'CA' DD' + (-) b 'CB f 

where x' is if X is a matrix of even forms or 1 if X is a matrix of odd forms. The 
associative property of this product can be easily checked. 
The supertrace is defined as 





StrM = Tr^ o ^JM Str I ^ ^ j = TrA - TtD (C.2) 

It is straightforward to prove the cyclic property of the supertrace, where the ^4-type 
(.F-type) supermatrices behave as odd (even) forms: 

Str(a b) = -Str(b a) ; Str(a f) = Str(f a) ; Str(f g) = Str(g f) . (C.3) 

We have denoted by o, b (f, g) generic supermatrices of the .A-type (F-type). 
We also define a pseudotransposition operation: 

(C.4) 

where t denotes usual matrix transposed. This generalizes the usual transposition in the 
sense that (using the same notation of (|C.3|) ): 




(a b) pt = -b pt a pt ; (a f) p * = f a pt ; (f Q ) pt = Q pt f pt ; (C.5) 

Nevertheless, unlike the usual transposition, pseudotransposing twice does not yield the 
initial matrix. Obviously, the supertrace does not change under pseudotransposition. 

D. Gauge transformation of the 1-form Q± 



In this appendix we show that the 1-form Q\ defined in ( |3.30| ) in section 3^ transforms non 



trivially only under U{\)a- This is an important point in the understanding of the U(1)a 
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axial anomaly and it is therefore worth being explicitly shown. For clarity of exposition, 
we report here the expression ( [3.30 ) for Cli 

= e - r2 Tr (iA L - iA R + (log T - log T^rdr) (D.l) 

It is obvious that the first two terms transform only under U(1)a transformations. The 
remaining log terms, instead, require a little more effort. First of all notice that we can 
write a general variation of Tr(logT — logT^) as 

5Tr(logT - logT 1 ") = Tr (t^ST - (T^StA = -^Tr (t ] 5T - T6tA = -^Tr (t*5T) 

(D.2) 

where in the second equality we used the same condition on T we imposed in section |3.3| , 
TT^ = T^T = T 2 Ij\f f , from where it also follows that T5T^ = —STT^, which we used in the 
last equality of (|D.2[) . 

From ( |A.7|) , we may show that for a vectorial (Aj, = A# = Ay) gauge transformation 
5\ v Tx(logT — logT^) vanishes. For an axial infinitesimal transformation (A^ = — Ar = 
Aa), from ( |A.7| ) and ( p.2| ), we find, instead: 

5 AA Tr(log T - log rt) = — 4TrA^ . (D.3) 

This is clearly non-zero for abelian ^(l)^ axial gauge transformations only, as we wanted 
to show. 



E. On the WKB approximation and linear confinement 



In order to use the WKB approximation to determine the spectrum, we start by map- 
ping the standard eigenvalue problem to a Schrbdinger-like equation. We follow pl| . We 
start from an equation (notice that A n is always proportional to the square of the four- 
dimensional mass): 



1 



T(z) 



Z 



+ B(z)i>{z) = X n ^(z) 



and define a new radial variable u and a rescaled wave function a: 



-r = S W 

az 



a(u) = E(u)ip(z(u)) , 



S(u) 



T(*(u)) 
E(*(u)) 



Then, equation (E.l) is rewritten as 12 : 

—a"(u) + V(u)a(u) = X n a(u) 

with: 

"(u)" 



(E.l) 



(E.2) 



(E.3) 



(E.4) 



12 In this appendix, we define V as the Schrodinger like potential that appears in equation (E.3). It should 
not be confused with the tachyon potential used in the rest of the paper 
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One can use the WKB approximation to estimate the mass of high excitations: 

-i 

(E.5) 



dX„ 

2vr 



dn 



1 du 

Ul 



^/X n -V(u)_ 

where u\ and u 2 are the classical turning points. In [^TJ it was observed that, quite 

generically, the range in the variable u is finite. Therefore, for large A n , the problem is 

_ i 

similar to an infinite potential well: the integral behaves as A n 2 and X n ~ to 2 ~ n 2 . In 



1 16] it was remarked that the behavior of the functions defined above can be tuned near 
the infrared such that the IR turning point in equation QE.5| ) goes to infinity in such a 
way that the desired relation m 2 ~ n is obtained. We may naively expect that highly 
excited masses should depend mainly on the UV behavior of different functions. However, 
linear confinement emerges from the IR dynamics and the size of highly excited hadrons 
grows with the excitation number. Therefore, highly excited mesons are still affected by 
IR physics. 



We now apply this formalism to the vector meson equation (4J3). Comparing to ( [E.lD , 



i 

we define X n = rr? n , B(z) = 0, T(z) = e~^gi z , S(z) = yff^r- The change of variable 
to u would involve a complicated integral, but in order to estimate the large excitation 
spectrum, it is enough to study the leading UV and IR behavior of the functions. 

Near the UV, we consider an asymptotically AdS space g Z z~ 9xx ~ ^AdSi ^ z ' > ^ ~ 
const, so S w 1 and ujjv = z - The potential (jE.4[) is Vuv( u ) ~ A- Thus, for large A n , 



the classical turning point is at u\ = -y/3/(4A n ) so u\ remains finite as X n — > oo and, as 

_i 

expected, the UV contribution to the integral in ( |E.5| ) decreases as X n 2 . 

The intermediate region when one cannot apply the IR nor the UV asymptotics has 

_ i 

finite size in the u- variable so its contribution to the integral also decreases as A n 2 . 

Near the IR, the vev of the open string tachyon is diverging and g xx Tqcd, 9zz ~ 
^(d z T) 2 and ~ e~ T . By substituting in ( |E.2| ), we find that the n-coordinate near 

the IR is just um ps \J ■kTqcd T ^ anc ^ ^ e P°t en ti a l is VlR ~ (f^QCD u) 2 . The classical 

i 

turning point is therefore at U2 = ^Tqcd ^" ' Thus, u 2 grows to infinity as A n , unlike the 

[is: 

+ . . . (E.6) 



cases considered in [41|. The IR contribution to the integral is: 

du 1 



where the dots stand for a piece that vanishes as A ra — » oo. Therefore, the integral in ( |E.5[ ) 
is dominated by the IR region and: 

lim ^ = 2ttT QC d (E.7) 

A n — >cxd an 



which, reinserting X n = yields ( 4.2Cj ) as we wanted to show. For the axial vectors, (see 



( 4.12 )), the z-dependent mass for A± (such that B(z) = g X xv(z) 2 ) adds an extra term to 
the potential near the IR Vir = (vr 2 /4 + A)T^ CD u 2 , leading to ( |4.21[) . 



Figure |2| depicts a sketch of the Schrodinger-like potential in the w-variable. 
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V(u) 




1 1 u 

Ui u 2 

Figure 2: Qualitative behavior of the Schrodinger-like potential in the u-variable. Near the UV, 
it grows as V oc u~ 2 while near the IR V oc u 2 . In the middle, it may present more complicated 
features which do not affect the leading behavior of the spectrum for large excitation number. 



References 



J. M. Maldacena, "The large N limit of superconformal field theories and super gravity," Adv. 
Theor. Math. Phys. 2, 231 (1998) [Int. J. Thcor. Phys. 38, 1113 (1999)] 



[ArXiv:hep-th/9711200] 



S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "Gauge theory correlators from 
non-critical string theory," Phys. Lett. B 428, 105 (1998) |[ArXiv:hep-th/9802109] 



[2] E. Witten, " Anti-de Sitter space and holography", Adv. Theor. Math. Phys. 2, 253 (1998); 
|[ArXiv:hep-th/9802150i . 

[3] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik and I. Kirsch, "Chiral symmetry 
breaking and pions in non- super symmetric gauge / gravity duals," Phys. Rev. D 69, 066007 
(2004) |[ArXiv:hcp-th/0306018l . 



[4] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, "Towards a holographic dual of 
[5] 



large-N(c) QCD," JHEP 0405 (2004) 041 |ArXiv:hep-th/0311270l 



N. J. Evans and J. P. Shock, "Chiral dynamics from AdS space," Phys. Rev. D 70, 046002 
(2004) |[ArXiv:hep-th/0403279] ; 



K. Ghoroku and M. Yahiro, "Chiral symmetry breaking driven by dilaton," Phys. Lett. B 



604, 235 (2004) |[ArXiv:hcp-th/0408040l 



N. Evans, J. Shock and T. Waterson, "D7 brane embeddings and chiral symmetry breaking," 



JHEP 0503, 005 (2005) [ArXiv:hcp-th/0502091] 



R. Apreda, J. Erdmenger and N. Evans, "Scalar effective potential for D7 brane probes which 



break chiral symmetry," JHEP 0605, 011 (2006) [ArXiv:hcp-th/0509219] 



[6] J. L. F. Barbon, C. Hoyos, D. Mateos and R. C. Myers, "The holographic life of the eta'," 
JHEP 0410, 029 (2004) |[ArXiv:hep-th/040426] . 

[7] T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD," Prog. Thcor. 



Phys. 113, 843 (2005) |[ArXiv:hep-th/04121411 



T. Sakai and S. Sugimoto, "More on a holographic dual of QCD," Prog. Theor. Phys. 114, 
1083 (2006) [ArXiv:hep-th/0507073l . 



-37- 



[9] 
[10] 

[11] 
[12] 
[13] 
[14] 
[15] 



O. Aharony, J. Sonnenschein and S. Yankielowicz, "A holographic model of deconfinement 



and chiral symmetry restoration, " [ArXiv:hcp-th/0604161] 



R. Casero, A. Paredes and J. Sonnenschein, "Fundamental matter, meson spectroscopy and 



non-critical string / gauge duality," JHEP 0601, 127 (2006) |[ArXiv:hep-th/0510110] 



E. Antonyan, J. A. Harvey, S. Jensen and D. Kutasov, "NJL and QCD from string theory, 



[ArXiv:hcp-th/0604017] 



K. Peeters, J. Sonnenschein and M. Zamaklar, "Holographic melting and related properties of 
mesons in a quark gluon plasma," Phys. Rev. D 74, 106008 (2006) |[ArXiv:hcp-th/0606195] . 



J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, "QCD and a holographic model of 



hadrons," Phys. Rev. Lett. 95, 261602 (2005) |[ArXiv:hep-ph/0501128] 



L. Da Rold and A. Pomarol, "Chiral symmetry breaking from five dimensional spaces," Nucl. 
Phys. B 721, 79 (2005) |[ArXiv:hcp-ph/0501218]| . 



J. P. Shock and F. Wu, "Three flavour QCD from the holographic principle," JHEP 0608, 



023 (2006) |[ArXiv:hcp-ph/0603142] 



C. Csaki and M. Reece, "Toward a systematic holographic QCD: A braneless approach, 



[ArXiv:hep-ph/0608266] 



N. Evans and A. Tedder, "Perfecting the ultra-violet of holographic descriptions of QCD, 



Phys. Lett. B 642, 546 (2006) [ArXiv:hcp-ph/0609 112] 



J. P. Shock, F. Wu, Y. L. Wu and Z. F. Xie, "AdS/QCD Phenomenological Models fr 



Back-Reacted Geometry, " | ArXiv :hcp-ph/0611227] 



[16] A. Karch, E. Katz, D. T. Son and M. A. Stephanov, "Linear confinement and AdS/QCD,' 



[ArXiv:hcp-ph/0602229] 



[17] 
[18] 

[19] 



U. Gursoy and E. Kiritsis, in preparation. 

S. Sugimoto and K. Takahashi, "QED and string theory," JHEP 0404, 051 (2004) 



[ArXiv:hcp-th/0403247] 



F. Bigazzi, R. Casero, A. L. Cotrone, E. Kiritsis and A. Paredes, "Non-critical holography 
and four- dimensional CFT's with fundamentals," JHEP 0510, 012 (2005) 



[ArXiv:hcp-th/0505140] 



[20] A. Sen, "Tachyon dynamics in open string theory," Int. J. Mod. Phys. A 20, 5513 (2005) 
|[ArXiv:hcp-th/0410103i . 

[21] Y. Kinar, E. Schreiber and J. Sonnenschein, "Q anti-Q potential from strings in curved 



spacetime: Classical results, " Nucl. Phys. B 566, 103 (2000) [ArXiv:hcp-th/9811192] 



J. Sonnenschein, "What does the string / gauge correspondence teach us about Wilson 



loops?," [ArXiv:hcp-th/0003032] 



[22] S. R. Coleman and E. Witten, "Chiral Symmetry Breakdown In Large N Chromodynamics, 
[23] 



Phys. Rev. Lett. 45, 100 (1980 



S. L. Adler, "Axial vector vertex in spinor electrodynamics," Phys. Rev. 177, 2426 (1969) 
J. S. Bell and R. Jackiw, "A PC AC puzzle: piO — > gamma gamma in the sigma model," 



Nuovo Cim. A 60, 47 (1969) 



[24] 



I. R. Klebanov and J. M. Maldacena, "Superconformal gauge theories and non-critical 
super strings," Int. J. Mod. Phys. A 19, 5003 (2004) |[ArXiv:hcp-th/0409133] . 



- 38 - 



[25] A. Fotopoulos, V. Niarchos and N. Prezas, "D-branes and SQCD in non-critical superstring 



theory," JHEP 0510, 081 (2005) |[ArXiv:hep-th/0504010l 



S. Murthy and J. Troost, "D-branes in non-critical superstrings and duality in N = 1 gauge 
theories with flavor," JHEP 0610, 019 (2006) |[ArXiv:hep-th/0606203]| 



A. Sen, "Dirac-Born-Infeld action on the tachyon kink and vortex," Phys. Rev. D 68, 066008 



(2003) [ArXiv:hep-th/0303057l 



M. R. Garousi, "D-brane anti-D-brane effective action and brane interaction in open string 



channel," JHEP 0501, 029 (2005) |[ArXiv:hep-th/0411222l 



E. Katz, A. Lewandowski and M. D. Schwartz, "Tensor mesons in AdS/QCD," Phys. Rev. D 
74, 086004 (2006) |[ArXiv:hep-ph/0510388]| . 



D. Kutasov, M. Marino and G. W. Moore, "Remarks on tachyon condensation in superstring 
field theory, " |[ArXiv:hep-th/0010108j . 



J. A. Minahan and B. Zwiebach, "Effective tachyon dynamics in superstring theory," JHEP 
0103, 038 (2001) [ArXiv:hcp-th/0009246l . 



A. A. Tseytlin, "On non-abelian generalisation of the Born-Infeld action in string theory," 
Nucl. Phys. B 501, 41 (1997) |[ArXiv:hcp-th/9701125l . 



M. R. Garousi, "Tachyon couplings on non-BPS D-branes and Dirac- Born-Infeld action, 
Nucl. Phys. B 584, 284 (2000) |[ArXiv:hep-th/0003122][ 



E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge 
theories," Adv. Theor. Math. Phys. 2, 505 (1998) |[ArXiv:hcp-th/98031311 . 



C. Vafa and E. Witten, "Restrictions On Symmetry Breaking In Vector - Like Gauge 



Theories, " |Nucl. Phys. B 234, 173 (1984) 



C. Kennedy and A. Wilkins, " Ramond-Ramond couplings on brane- antibrane systems" , 
Phys. Lett. B 464, 206 (1999); |[ArXiv:hcp-th/9905195] . 



P. Kraus and F. Larsen, " Boundary string field theory of the DD-bar system" , Phys. Rev. D 
63, 106004 (2001); |ArXiv:hep-th/0012198l . 



T. Takayanagi, S. Terashima and T. Uesugi, " Brane- antibrane action from boundary string 



field theory", JHEP 0103, 019 (2001); |[ArXiv:hcp-th/0012210l 



M. B. Green, J. A. Harvey and G. W. Moore, "I-brane inflow and anomalous couplings on 
D-branes," Class. Quant. Grav. 14, 47 (1997) |[ArXiv:hcp-th/9605033]. 



D. Z. Freedman, S. D. Mathur, A. Matusis and L. Rastelli, "Correlation functions in the 
CFT(d)/AdS(d + 1 ) correspondence, " Nucl. Phys. B 546, 96 (1999) |[ArXiv:hep-th/9804058][ 



R. Jackiw, "Topological Investigations of Quantized Gauge Theories, " in Current Algebra and 
Anomalies, eds. S. B. Treiman, R. Jackiw, B. Zumino, and E. Witten, ed. World Scientific 
Publishing, 1985. 



E. Schreiber, "Excited mesons and quantization of string endpoints," [ArXiv:hep-th/0403226] 



M. Shifman, "Highly excited hadrons in QCD and beyond," [ArXiv:hcp-ph/0507246]| . 



M. Gell-Mann, R. J. Oakes and B. Renncr, "Behavior of current divergences under SU(3) x 
SU(3)," Phys. Rev. 175, 2195 (1968). 



- 39 - 



[44] A. Armoni, "Witt en- Ven eziano from Green- Schwarz," JHEP 0406, 019 (2004) 
|[ArXiv:hcp-th/0404248i . 

[45] O. Bergman and G. Lifschytz, "Holographic U(1)A and string creation," 



[ArXiv:hep-th/0612289] 



J. M. Maldacena and C. Nunez, "Towards the large N limit of pure N = 1 super Yang Mills," 



Phys. Rev. Lett. 86, 588 (2001) [ArXiv:hcp-th/0008001] 



[47] 

[48] E. Witten, "Current Algebra Theorems For The U(l ) Goldstone Boson, " iNuci. Phys. B 156 



E. Witten, "Theta dependence in the large N limit of four- dimensional gauge theories," Phys. 
Rev. Lett. 81, 2862 (1998) |[ArXiv:hep-th/9807109]| . 



269 (19791 ; 



G. Veneziano, "U(l) Without Instantons," fhicl. Phys. B 159, 213 (1979J . 

L. Da Rold and A. Pomarol, "The scalar and pseudoscalar sector in a five-dimensional 



approach to chiral symmetry breaking," JHEP 0601, 157 (2006) [ArXiv:hcp-ph/0510268] 



S. Kuperstein and J. Sonncnschcin, "Non-critical supergravity (d > 1) and holography,' 



JHEP 0407, 049 (2004) [ArXiv:hep-th/0403254] 



S. Kuperstein and J. Sonncnschcin, "N on- critical, near extremal AdS(6) background as a 
holographic laboratory of four dimensional YM theory," JHEP 0411, 026 (2004) 



[ArXiv:hcp-th/0411009] 



[51] D. Mateos, R. C. Myers and R. M. Thomson, "Thermodynamics of the brane," 
|[ArXiv:hep-th/07011~32] . 

[52] R. Casero, C. Nunez and A. Paredes, "Towards the string dual of N = 1 SQCD-like 



theories," Phys. Rev. D 73, 086005 (2006) |[ArXiv:hcp-th/0602027]| ; 

A. Paredes, "On unquenched N = 2 holographic flavor," JHEP 0612, 032 (2006) 



[ArXiv:hcp-th/0610270] 



F. Benini, F. Canoura, S. Cremonesi, C. Nunez and A. V. Ramallo, "Unquenched flavors in 



the Klebanov- Witten model," [ArXiv:hcp-th/0612118] 



[53] K. Skenderis, "Lecture notes on holographic renormalization," Class. Quant. Grav. 19, 5849 
[54] 



(2002) |[ArXiv:hep-th/0209067] 



A. Karch, A. O'Bannon and K. Skenderis, Holographic renormalization of probe D-branes in 
AdS/CFT, " JHEP 0604, 015 (2006) |[ArXiv:hep-th/0512125] . 



- 40 - 



